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We develop an analytical approach for studying lattice gauge theories 
within the plaquette representation where the plaquette matrices play the 
role of the fundamental degrees of freedom. We start from the original Ba- 
■ trouni formulation and show how it can be modified in such a way that each 

oo : non-abelian Bianchi identity contains only two connectors instead of four. In 

addition, we include dynamical fermions in the plaquette formulation. Using 
this representation we construct the low-temperature perturbative expansion 
for U(l) and SU(N) models and discuss its uniformity in the volume. The 
final aim of this study is to give a mathematical background for working with 
O ■ non-abelian models in the plaquette formulation. 



1 Introduction 



- 1—1 

X. 

1.1 Motivation 



There exist several equivalent representations of lattice gauge theories (LGT). Orig- 
inally, LGT was formulated by K. Wilson in terms of group valued matrices on links 
as fundamental degrees of freedom Xj. The partition function can be written as 

DU exp{(3S[U n (x)}} , (1) 

where S[U n (x)} is some gauge-invariant action whose naive continuum limit coincides 
with the Yang-Mills action. The integral in (0) is calculated over the Haar measure 
on the group at every link of the lattice. Very popular in the context of abelian 
LGT is the dual representation which was constructed in [2]-|I|- Extensions of 
dual formulations to non-abelian groups have been proposed only in the nineties in 
[5]-[TU]. The resulting dual representation appears to be a local theory of discrete 
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variables which label the irreducible representations of the underlying gauge group 
and can be written solely in terms of group invariant objects like the 6j-symbols, 
etc. A closely related approach to the dual formulation is the so-called plaquette 
representation invented originally in the continuum theory by M. Halpern [TT] and 
extended to lattice models in [T2|. In this representation the plaquette matrices play 
the role of the dynamical degrees of freedom and satisfy certain constraints expressed 
through Bianchi identities in every cube of the lattice. Each representation has its 
own advantages and deficiencies. E.g., the Wilson formulation is well suited for 
Monte-Carlo simulations, while dual and plaquette representations are usually used 
for an analytical study of the models. In particular, duals of abelian U(l) LGT 
have been used to prove the existence of the deconfinement phase transition at 
zero temperature in four- dimensions (4D) fSl an d to prove confinement at all 
couplings in 3D Also Monte-Carlo simulations proved to be very efficient in 
the dual of AD U(l) LGT fl5|. 

So far, however both dual and plaquette formulations have not been so pop- 
ular in the case of non-abelian models, probably due to the complexity of these 
representations. For instance, the plaquette representation can hardly be used for 
Monte-Carlo computations due to a number of constraints on the plaquette ma- 
trices. Let us remind the general form of the Batrouni construction . In ^2] the 
plaquette representation was constructed in the maximal axial gauge. The partition 
function takes the following form if S[U n (x)] in (JTJ is the standard Wilson action 



where V p G SU(N) are plaquette matrices, DV P is the invariant Haar measure of 
SU (N) . The product over c runs over all cubes of the lattice. The Jacobian J(V C ) 
is given by 



where the sum over r is a sum over all representations of SU(N), d r = Xr{I) is the 
dimension of the representation r. The last expression is nothing but an SU(N) 
delta-function which introduces certain constraint on the plaquette matrices. This 
constraint is just the lattice form of the Bianchi identity. The SU(N) character 
Xr depends on an ordered product of SU(N) matrices as dictated by the Bianchi 
identity. Its exact form will be given in the next section. An important point is that 
in the non-abelian case the resulting constraints on the plaquette matrices appear 
to be highly-nonlocal and this fact makes an analytical study of the model rather 
difficult. In particular, it has prevented so far the construction of any well controlled 
and useful weak-coupling expansion. 

A different plaquette formulation of 3D SU(N) LGT has been proposed in [T7| . 
It has a local form and does not require gauge fixing. Unfortunately, as we have 
found by explicit computations the model proposed in J7| does not coincide with 
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the Wilson LGT contrary to the claim of ^7j. Exact calculations on a 3 3 lattice 
show the equivalence between Wilson LGT and the model of [T7] but the equivalence 
is lost already for a 4 x 3 2 lattice. The point is that the constraints on the plaquette 
matrices in the model of do not match the non-abelian Bianchi identity as will 
be seen from our explicit calculations. Nevertheless, we have found that a certain 
decomposition of the lattice made in ^7J can be useful in simplifying Batrouni's 
original representation. In the next section we use some ideas of ^7J to reduce 
the number of connectors in constraints on plaquette matrices from four to two per 
each cube of the lattice. This simplifies the whole representation but still it is quite 
involved. 

Let us also mention that there exists a plaquette representation in terms of so- 
called gauge-invariant plaquettes [18]. This representation does not require gauge 
fixing and can be formulated both on finite and infinite lattices. It is quite possible 
to work also with this formulation, all the methods developed in this paper can 
be straightforwadly extended to the model of gauge-invariant plaquettes. We have 
nevertheless found that the plaquette formulation obtained in the maximal axial 
gauge is simpler to handle, especially on finite lattices. Moreover, we shall explain 
how our formulation can be extended to periodic lattices. In addition to the previous 
works |T2|, JJJ, JEj we include also dynamical fermions in the plaquette formulation. 

In spite of the complexity of the plaquette representation, we think it has certain 
advantages compared to the standard Wilson representation. Some of them have 
been mentioned and elaborated in ^2] and ^H] • Duality transformations, Coulomb- 
gas representation, strong coupling expansion look more natural and simpler in the 
plaquette formulation. It is also possible to develop a mean-field method which is 
gauge-invariant by construction and is in better agreement with Monte-Carlo data 
than any mean- field approach based on the mean-link method [T§] . 

Nevertheless we believe that the main advantage of this formulation, not men- 
tioned in ^2] 5 |JE] and ^H] lies in its applications to the low-temperature region. 
Let V p be a plaquette matrix in SU (N) LGT. The rigorous result of [[SO] asserts that 
the probability that Tr(J — V p ) > £ is bounded by 

p(0 < 0(e~ m ) , $ -»• oo , b = const (4) 

uniformly in the volume. Thus, all configurations with £ > 0(f3~ l ) are exponentially 
suppressed. This is equivalent to the statement that the Gibbs measure of SU(N) 
LGT at large (3 is strongly concentrated around configurations on which V p ~ I. 
This property justifies expansion of the plaquette matrices around unity when (3 is 
sufficiently large while there is no such justification for the expansion of link matrices, 
especially in the large volume limit. In particular, we think that replacing TrV p in 
the Gibbs measure by a Gaussian distribution, in the region of sufficiently large f3 
is a well justified approximation. In fact, all the corrections to this approximation 
must be non-universal. 

Actually, this is one of our motivations to construct a low-temperature, i.e. large- 
f3 expansion of gauge models using the plaquette representation. The well-known 
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problem of the standard perturbation theory (PT), i.e. whether PT is uniformly 
valid in the volume, can be shortly formulated in the following way. When the 
volume is fixed, and in the maximal axial gauge the link matrices perform small 
fluctuations around the unit matrix and the PT works very well producing the 
asymptotic expansion in inverse powers of f3. However, in the large volume limit the 
integrand becomes arbitrarily flat even in the maximal axial gauge. It means that in 
the thermodynamic limit (TL) the system deviates arbitrarily far from the ordered 
perturbative state, so that no saddle point exists anymore, i.e. configurations of 
link matrices are distributed uniformly in the group space. That there are problems 
with the conventional PT was shown in [2T] . where it was demonstrated that the PT 
results depend on the boundary conditions (BC) used to reach the TL. Fortunately, 
even in the TL the plaquette matrices are close to unity, the inequality (jlj) holds 
and thus provides a basis for the construction of the low-temperature expansion 
in a different and mathematically reliable way. In this paper we develop such a 
weak-coupling expansion both for abelian and non-abelian models. 

This paper is organised as follows. In the next section we give our plaquette for- 
mulation of 3D SU(N) LGT. We work in maximal axial gauge and consider a model 
with arbitrary local pure gauge action. For fermions we choose either the Wilson 
or the Kogut-Susskind action. The plaquette representation will be formulated on 
a dual lattice for the partition function, 't Hooft and Wilson loops. In section 3 we 
construct the weak-coupling expansion for the abelian model using the plaquette 
representation. We give a general expansion for the partition function, calculate the 
zero-order generating functional and show how to compute corrections and expec- 
tation values of Wilson loops. Then, we extend the weak-coupling expansion to an 
arbitrary SU(N) gauge model. Here we give a general expansion of the Boltzmann 
factor, explain how to treat the Bianchi constraints in the expansion, compute the 
generating functional and establish some simple Feynmann rules. Finally, we dis- 
cuss some features of the large-/? expansion in non-abelian models. Our conclusions 
are presented in section 4. Some computations are moved to the Appendices. In 
Appendix A we study the link Green functions which appear as the main building 
blocks of the expansion in the plaquette formulation. In the Appendices B and C we 
give all technical details for the calculation of the free energy expansion for SU(N) 
models in the plaquette representation. 

1.2 Notations and conventions 

We work on a 3D cubic lattice A G Z d with lattice spacing a = 1, a linear extension 
L and x = (x,y,z), x,y,z G [0, L — 1] denote the sites of the lattice. We impose 
either free or Dirichlet BC in the third direction and the periodic BC in other 
directions. Let G = U(N),SU(N); U t = U n (x) G G, V p G G and DU h DV P denote 
the Haar measure on G. Xr{U) and d(r) will denote the character and dimension 
of the irreducible representation {r} of G, correspondingly. We treat models with 
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local interaction. Let H[U] be a real invariant function on G such that 

I Hp] I < H{I) (5) 
for all U and coefficients of the character expansion of exp(/3H) 

C[r] = J DU exp((3H[U]) X r(U) (6) 

exist. Introduce the plaquette matrix as 

V p = U n (x)U m (x + e n )Ul(x + e m )Ul{x) , (7) 

where e n is a unit vector in the direction n. The action of the pure gauge theory is 
taken as 

S g [U n {x)] = Y,H[V p ] . (8) 

peA 

The action for fermions we write down in the form (colour and spinor indices are 
suppressed) 

N f 

x,x'£A f=l 

where 

1 d 

A f (x, x'; Ut) = M f 5 X;X , + [5 x+en , x 'Ux)Un(x) + 5 x - en , x ,Z n (x')Ul(x')] . (10) 

A n=l 

We have introduced here the following notations 

M f = m f -rd , £ n (x) = r + 7„ , ^ n (x) = r - 7„ (11) 

for Wilson fermions and 

M f = m f , e„(x) = = ^n(^) = (12) 

for Kogut-Susskind fermions. m/ is mass of the fermion field, iV/ is the number of 
quark flavours and r is the Wilson parameter (r = 1 is a conventional choice). 

After integrating out fermion degrees of freedom the partition function of the 
gauge theory on A with the symmetry group G can be written as 

Z A (/3,m f ,N f )= /n^^xexpL^Z/^+^Trln^^x';^)! , (13) 

x,n [ p G A f=l J 

where the trace is taken over space, colour and spinor indices. 
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2 Plaquette formulation and expression for the 
Jacobian 



In this section we give our formulation of the plaquette representation for 3D SU(N) 
LGT. A short description of our procedure can be found in [22] for pure gauge theory 
and in (22] for a theory with fermions. In subsection 2.1 we calculate the plaquette 
representation for the partition function on a dual lattice. In the second subsection 
we derive the plaquette representation for some observables. 

2.1 Expression for the Jacobian 

Let us first give a qualitative explanation of our transformations. Consider the 
partition function given in Eqs. (j2J) and @. For the formulation of the Bianchi 
identity related to a given 3(i-cube we define according to 17J a vertex A in this 
cube and separated by the diagonal a vertex B, see FigUJ This assignment we extend 
to all neighbouring cubes and finally to the whole lattice. Thus, all A vertices are 
separated by 2 lattice spacings, the same is true for B vertices. Next, we take a 
path connecting the vertices A and B by three links U±, U 2 , U3 as shown in FigO] 
Then the matrix V c in (J2J) and entering the Bianchi identity for the cube c can 
be presented in the following form 



FJ A means an appropriately ordered product over three plaquettes of the cube 
attached to the vertex A. The matrix C defines a parallel transport of vertex A 
into vertex B and equals the product of three link matrices connecting A and B, 
see FigJT] 



The connector C plays a crucial role in the non-abelian Bianchi identity. Its path has 
to fit to the choices of YI p <=a an< ^ Ftpes- ^ e cnoose t ne structure of the connectors 
as shown in Fig|2] This 8-cube fragment of the lattice is repeated through the whole 
lattice by simple translation. As is seen from Fig|2] there are only four different 
types of connectors. For example, the connector B\A is the same as B 7 A. The 
collection of cubes with the same type of connectors, e.g. B\A form on a dual 
lattice a body centred cubic (BCC) lattice with double lattice spacing. There are 
four different sub-lattices of this type corresponding to the four types of connectors 
in Fig|2] B±A, B2A, B3A and B4A, and therefore four different types of Bianchi 
identities. 

Consider now the partition function given by Eq. (|13|) on a 3D lattice with free 
or Dirichlet BC in the third direction and periodic BC in other directions for the 
link matrices. To get the plaquette representation we make a change of variables 




(14) 



c = u x u 2 ul . 



(15) 
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Figure 1: Vertices A, B and the connector path on the cube. 




Figure 2: Structure of connectors on the lattice. 
(JZJ) in the partition function fH?|) . Then the partition function gets the form 

Z A (/3,m f ,N f ) = j n^expj^^F^m J(V P ) , (16) 



where the Jacobian of the transformation reads 



J(V P ) = f ftdUt ]J6 ( V; ^ - l\exp I ^Tr In A 
J i p \ lep J I f=i 



f (x,x'; Ui) 



(17) 



The last equation is rather formal, for we have to specify the order of multiplications 
of non-abelian matrices. An important point concerns the position of the plaquette 
matrix within the product of link matrices. The plaquette matrices V p we insert 
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at the vertices A or B of each plaquette attached to A or B, correspondingly. To 
write down the plaquette delta-functions in (|17|) explicitly we use the following 
conventions: 1) a plaquette is always defined as U n (x)U m (x + e n )U^(x + e m )L r J ^(x) 
with n < m, i.e. it starts from the smallest possible coordinate; 2) plaquette matrices 
are inserted as V p , not VJ in all vertices A and B. Then, for instance, the six delta- 
functions for the cube B\A in FigOcan be written with the notations of Fig|3]as 

6 (v pi u h u l2 ulul - 1) , s(y P2 u h u h ulul-i) , 

5 (v m U k U l7 UlX -i) ,5 (u l3 U l6 V P4 Ulul - /) , (18) 

6 (u h U k V P5 Ul o Ul -i) ,8 (u l9 U h0 V P6 Ulul - i) . 

The change of variables in (J7J) is uniquely determined from these expressions. Similar 
expressions can be written for all types of connectors using the same rules. Now we 
want to map one of the plaquette delta's to a delta-function for the cube using all 
other delta's. This procedure can be done for all cubes of the lattice in the following 
way. We start from the plane z = L — 1 and use space-like plaquettes (i.e., lying 
in the xy-pl&ne) to map them onto cube delta-function. For the above cube B\A 
it means that we take the plaquette delta for V P6 and substitute link matrices by 
their expressions obtained from all other delta's. Then we repeat these calculations 
for all cubes. One should be careful about the order in which the cubes are taken 
since to accomplish such calculations for all cubes, the order is important. We can 
proceed strictly plane by plane and stop at the plane z = 1. There are still space- 
like plaquettes left in the plane z = 0. The connector for a given cube is taken in 
accordance with Fig|21 Therefore, we mapped all space-like plaquette delta's lying 
on the top of a given cube onto cube deltas with connectors lying on the bottom of 
this cube (in the gauge U 3 = I defined below). For instance, if the original plaquette 
lies in the plane z = Zi, then the connector lies in the plane Zi — 1. These cube delta- 
functions become samples for the Bianchi constraints. In this way we arrive at the 
following expression for the Jacobian (|17j) 

j(v P ) = f n du i n 6 ( v po n ^ - j ) n n j{ y^ 

l po \ iepo / »=l c(i) 

( Nf ) 

xexp < ^TrlnAf&x';^) > , (19) 

where J(V C ) is given in (jSJ). n po runs over a ^ time-like plaquettes and Ylctt) means 
a product over all cubes with the i-th type of connector. In the unique notations of 
FigEH the matrices Vc for different types of connectors Cb a (see FigJ2J) are given 
by 

Vc = V p % 2 V p \C BlA V P5 V P6 V p \C BiA , C BlA = U h U h U k , (20) 
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v c 2 = v p \v p [v p \c B2A v P2 v P6 v P3 c B2A , c B2A = u h u}u k 



(21) 



v c = v P5 v p \v p \c B3A v p \v P6 v P2 c B3A , c BzA = u{u}u h , 
Vc = v P2 v P5 v p \c BiA v p {v P6 v p {c B4A , c BiA = ulu u u h . 
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(22) 
(23) 



Figure 3: Notations for links and plaquettes of the original cube B\A. 

The rest of the derivation follows closely Ref.p2]- Choose the maximal axial 
$auge in the form 



U 3 {x,y,z) = U 2 {x,y,0) = U 1 {x,0,0) 



(24) 



For the Dirichlet BC one has ^(x, y, 0) = Ui(x, y, 0) = /. Note, that the gauge fix- 
ing procedure and the procedure described above are freely interchangeable. After 
gauge fixing one can easily obtain expressions for link variables in terms of pla- 
quettes. First of all, space-like links from the plane z = L — 1 can be trivially 
integrated out removing all plaquette but Bianchi delta's from the corresponding 
cubes near the top boundary of the lattice. Further, we observe that expressions 
for link variables will depend on the position of a link relatively to the vertex A or 
B. As example, consider a line containing vertices B in z direction, x, y being fixed 
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and let one of the vertices B have the coordinates B = (x, y, z). Then one finds for 
the link I = (x, y, z; n 2 ) 

o 

U n2 (x,y,z)= Yl Vmfay,^) , (25) 

z'=z-l 

while for the link I — (x,y — e 2 , z\ n 2 ) 

z-l 

U m (x,y-e 2 ,z) = JJ V 23 (x,y - e 2 ,z') , (26) 

z'=0 

where V23 is the plaquette matrix from the (zy)-plane. Similar expressions one 
obtains for all links entering in or coming out of the considered line. Obviously, the 
same rules apply for links in direction n\ as well as for lines containing vertices A. 
Using these expressions one can easily integrate out all link variables in (jl9j) . Clearly, 
the presence of the fermion determinant cannot change the integration procedure. 
In what follows we pass to the dual lattice where cubes become sites, links become 
plaquettes and sites become cubes. On the dual lattice our original notations of Fig|3] 
for the cube B\A become as in Fig0] This helps us to obtain explicit expressions 
for the Jacobian. 

Now we can write down the plaquette representation for the SU(N) gauge model 
on the dual lattice 

I f=l 

4 

x n n j w°) > ( 2? ) 

i=l x(i) 

where i = 1,2,3,4 denote four BCC sub-lattices, Yl x (i) runs over au s ^ es °^ sucn 
a BCC sub- lattice and J(Vx^) is an SU(N) delta- function (J3J). To visualise calcu- 
lations in the next sections we give here full expressions for the matrices ■ As 
example, we also give the graphical representation of the Bianchi identity for the 
cube BiA in Fig El To write down the full expressions we choose a coordinate sys- 
tem on the dual lattice as shown in FigEl Let afj = (xi, y^ Zi) be a site of the dual 
lattice which belongs to the i-th BCC sub-lattice and introduce notations for dual 
links 

h = (£J;ni) , l 2 = (xi;n 2 ) , l 3 = (^;n 3 ) , 
h = (xi - ex] rii) , l 5 = (xi - e 2 ; n 2 ) , l 6 = (xl - e 3 ; n 3 ) . (28) 

We keep these notations for links attached to any site x^. Then, for the type of the 
connector for the cube B\A one has 

= vlv h vl c m v l2 v l3 vl d* , (29) 



Z K {p,m f ,N f ) = 
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Figure 4: Original cube B\A on the dual lattice. 

1 Xi 2j — 1 

^(i) = II v n 2 {xu Vi ~ e 2 , JJ K 3 (j, 0) JJ V ni (xi - ei, , 

k=Zi — 1 j=L—l p=l 

for the connector of the cube -B 2 ^4 




1 L-l 2,-1 

k=Zi — 1 j=%i P=l 

for the connector of the cube 5 3 A 




1 Xj + l Zj-l 

^(3)= II W»*0 II Vl(j : y i: 0)l[vl(x i: y i:P ) , 

k=Zi-l j=L—l p=l 
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X 



Figure 5: Graphical representation on the dual lattice of the cube delta-function 
for the cube B\A. 

and for the connector of the cube B^A 

= V h V l2 Vl C m VlV h Vl C| (4) , (35) 

1 L-l Zi-1 

(xi, Vi - e 2 ,p) • (36) 

k=Zi-\ j=Xi+l p=l 

Finally, the fermionic matrix in the gauge (J2"2j) can be presented on the dual lattice 

as 

A f (c(x), c(x'); UtiVi]) = M f 6 c>c , + - [5 c+e3 , c ,£ 3 (c) + 5 c „ e3 ^ 3 (c')] 
1 2 

+ 2 E [<We„, c 'ac)^[c; ^] + ^^(c^c'; ^]] , (37) 

n=l 

where c(x) is a cube dual to the site x. U n [c; Vj\ is given in (|2"5j) . (j2Ej) for n = 2 and 
similar expressions can be easily written for n — 1. 

The expressions for the connectors C^(j) in (|29j) -(|3fi j) are given for the free BC. 
For the Dirichlet BC one should omit the second product ]~{ . V n3 (j) from formulae 
(jSOjl . (J32*|) . (J53j) and (J5fi|) . Extension of this representation to the lattice with periodic 
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BC in all directions is, at least in principle very simple. Here we give a qualitative 
explanation only, the details will be given elsewhere. Instead of the maximal axial 
gauge which is not compatible with the periodic BC one can fix the following gauge 

L-l 

U 3 (x, y, z) = / , z e [1, L - 1] ; U 3 (x, y, 0) = W(x) = J] U 3 (x, y, z) , (38) 

2=0 

where W is the Polyakov line. All the steps performed above are unchanged except 
for the integration over space-like links in the time slice z — 0. In this slice the 
expressions for connectors are modified: they include the values of the Polyakov 
lines. Therefore, one can express the final result in terms of plaquettes and these 
Polyakov lines (|38j). 



2.2 Observables on the dual lattice 

In this subsection we derive on the dual lattice a representation for the t' Hooft line 
and the Wilson loop. 

We consider first a t' Hooft line which in 3D consists of a string S xy of dual links 
between the sites x and y on the dual lattice. For simplicity let us restrict ourselves 
to pure gauge theory with the Wilson action. With the set T of plaquettes dual to 
the string S xy we define the action 

S Z (U; T)=(3J2 R eTr ^P + P ReTrzU P ( 39 ) 

and the partition function 

Z(z; S xy ) = j DU exp{(3S z (U;T)} , (40) 
where z is center element of SU(N). We can introduce the disorder operator as 



= Wq ■ (41) 



On the dual lattice the partition function Z(z, S xy ) of pure gauge theory can be 
written in the form of Eq. (j27|) with the action 

S Z (V; S xy ) ={3J2 ReTrV ^ + 1 3 Yl ReTlzV i ■ ( 42 ) 

Suppose that the points x and y belong to the BCC sub-lattice with the i — 1 type 
of connector. By a change of variables VJ — > z*Vi for I 6 S xy we get the following 
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representation for the partition function Z(z\ S, 



xyj 



Z(z; S xy ) = J \\dVi exp 



nn ^w°) 

i=l x(i) 



n hvp) 



J(zVP)J(zVU) 



(43) 



since z cancels from all Bianchi constraints but those for the endpoints. Therefore, 
one can say that the Bianchi constraint is violated at these endpoints. For example, 
for the SU(2) gauge group z = — 1 and J(zV) has the form 



j{zV) = 8{V + !) = J2(-l) 2r d rXr (V) 



(44) 



This delta-function constrains the product of non-abelian matrices to —I. 

Now we give a dual expression for a Wilson loop of size R x T in some repre- 
sentation j. Suppose, for simplicity that the loop contour lies in the y — z plane 
and R, T are even. Generelizing the consideration of the previous subsection it is 
straightforward to show that such a Wilson loop has the following dual form 

R/2-1 / z-1 

W j (C) = (Ti j J] J] V 1 (x,y + 2n,z 1 )l[V 1 (x,y + 2n + l,z 2 ) 

n=0 \zi=z— 1 Z2=0 
/z+T-l \ 

II II VHx,y + 2n + l, Zl ) l[ V?(x,y + 2n,z 2 ) ). (45) 

n=R/2-l V zi=0 z 2 =z+T-l / 

As one can see, the connectors are not cancelled in the general from the Wilson 
loop. One can simplify this expression taking one side of the loop in the plane 
z = 0. Then the Wilson loop is an ordered product of dual links which all lay inside 
the area spanned by the loop C, i.e. the first line in the last formula is unity. 



3 Low-temperature expansion 

This section is devoted to the construction of the low-temperature expansion of 
gauge models in the plaquette formulation. In doing this we follow the strategy 
developed for the two-dimensional principal chiral models in [21]. In that paper 
a low-temperature expansion was constructed for SU(N) spin models in the link 
formulation. This formulation appears to be closely related to the plaquette formu- 
lation of 3D gauge models. For example, as was shown in ^2] and i n both cases 
the strong coupling expansion is the expansion towards restoration of the Bianchi 
identity. As will be seen below, the weak-coupling expansion in both models also 
bears many common features. In the abelian case they turn out to be practically 
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identical, while the only (but essential) difference in the non-abelian gauge models 
is the appearance of connectors. 

A lattice PT in the maximal axial gauge was constructed for abelian and non- 
abelian models in the standard Wilson formulation by F. Miiller and W. Riihl in [26 . 
In this gauge the PT for non-abelian models displays serious infrared divergences in 
separate terms of the expansion starting from 0{f3~ 2 ) order for the Wilson loop and 
from 0(f3~ l ) for the free energy. Of course, it is expected that all divergences must 
cancel in all gauge invariant quantities. In [2B] the authors worked out a special 
procedure of how to deal with such divergences. The essential ingredients of the 
procedure of are the following: 1) fixing the Dirichlet BC in one direction and 
the periodic ones in other directions; 2) a Wilson loop under consideration is placed 
at distance R from the boundary and the limit R — > oo is imposed to restore time 
translation invariance; 3) divergent Green functions must be properly regularised, 
though there is no a priori preference to any regularization. With this procedure the 
first two coefficients of the expansion of Wilson loop expectation values in £77(2) 
LGT have been computed. These coefficients appear to be infrared finite in all D > 3 
and are proportional to the perimeter of the loop in 4D and to R In T + T In R in 
3D (for a rectangular loop R x T) . 

We are not aware of any proof of the infrared finiteness at higher orders of the 
expansion. That there could be a problem with the infrared behaviour was shown 
however in [2*T] . Namely, in the addition to the Dirichlet BC on the boundary, one 
link was fixed in the middle of the lattice and it was shown that the TL value of the 
expectation value of the plaquette to which the fixed link belong differs from that 
obtained in The underlying reason of such infrared behaviour lies, of course in 
the fact that the conventional PT is done around the vacuum U n (x) = I which is 
the true ground state only at fixed L, even in the maximal axial gauge. When the 
volume grows the fluctuations of the link matrices become larger and larger and may 
cause an infrared unstable behaviour. In particular, the integration regions over fluc- 
tuations are usually extended to infinity. But only for fixed L one can really prove 
that this introduces exponentially small corrections. There is no tool for proving 
that they remain exponentially small also in the large volume limit. On the contrary, 
bounds on the large plaquette fluctuations (J3J holds uniformly in the volume imply- 
ing that large plaquette fluctuations remain exponentially small also in the TL. This 
is the first achievement of the low-temperature expansion in the plaquette formula- 
tion. Constructing the low-temperature expansion in the plaquette formulation we 
aim not only to develop a new technics of calculation of the asymptotic expansions 
in LGT but also to get a deeper insight into this infrared problem and as a conse- 
quence into the problem of the uniformity of the low-temperature expansion. We 
hope that this will help in the investigation of the general properties of asymptotic 
expansions of non-abelian gauge models. As will be seen below, all Green functions 
appearing in the plaquette formulation are infrared finite. The source of trouble are 
the connectors of the non-abelian Bianchi identity. The low-temperature expansion 
in the plaquette formulation starts from the abelian Bianchi identity (zero order 
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partition function) for every cube and goes towards gradual restoration of the full 
identity with every order of the expansion. Thus, the generating functional contains 
an abelianized form of the identity without connectors. High-order terms include 
expansion of connectors (of course, together with other contributions which are how- 
ever infrared finite) which can lead to the infrared problem. Let us also mention that 
in this respect there is a big difference between gauge models and 2D SU(N) spin 
models. In the latter the low-temperature expansion also goes towards restoration of 
the full non-abelian Bianchi identity [21] and the form of the generating functional 
is formally the same. It includes invariant link Green functions and the standard 
Green function. In 2D the latter is not infrared finite and is the main source of the 
trouble. As will be seen below, in 3D gauge models all Green functions entering 
the generating functional are infrared finite. The source of the infrared problem in 
gauge models are the connectors of the non-abelian Bianchi identities. 

3.1 U(l) LGT 

For a simple introduction to our method we first consider the abelian U(l) LGT 
without fermions where the expansion can be done in a straightforward manner. 
Due to a cancellation of the connectors from Bianchi identities for abelian models 
the plaquette formulation for the U(l) model on the dual lattice with the Wilson 
action reads 

n>> ( 46 ) 

X 

where the Jacobian is given by the periodic delta-function 

oo 

Jx = Yl eir " x ' ( 4? ) 

r=— oo 

u x = uj ni (x) + u n2 (x) + uj n3 (x) - u ni (x - ei) - u n2 (x - e 2 ) - uJ n3 (x - e 3 ) . 

Dual links are defined by the point x and the positive direction n, see Eq. ()28jl . 
Strictly speaking, if the original link gauge matrices satisfy free boundary conditions, 
and so do the plaquette matrices then the dual links ui and representations r x must 
obey zero Dirichlet BC. All general expansions given below are valid for any type 
of BC, the dependence on the BC enters only the Green functions defined below. 
In 3D the difference between Green functions for different types of BC is of the 
order 0(L~ 2 ) for large L. Since we are interested in the TL behaviour only, we can 
consider dual models without any reference to the original link representation and 
introduce any type of BC. Both in abelian and in non-abelian cases we work with 
periodic BC. We remind at this point that expressions obtained below on a finite 
lattice do not correspond to the standard link formulation since the latter would 
include contributions from nontrivial Polyakov loops on a periodic finite lattice. 
However, in the TL both models must coincide and convergence to the TL is very 
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fast, like 0(L~ 2 ). If one wants to recover the exact correspondence between our 
dual representation and the link representation on the finite lattice with free BC, 
the only change is that one should take the corresponding Green functions entering 
the generating functionals. 

The construction of the PT proceeds as follows. The first step is a standard one, 
i.e. we re-scale u —>■ uj \/j3 and expand the Boltzmann factor in powers of 1/(3 as 



cxp 



(3 cos 



VP. 



exp 



l+J2((3)- k A k (u;f) 



k=l 



(4£ 



where are known coefficients (see, e.g. [23 )• Introducing sources hi for the link 
field, we find 



Z u{1) ((3»l) = eW L3 - L3 ^H 



OO 



fc=l 



I (I 2 

(3 k k [dh 2 



M u{1) (hi) , 



(49) 



where Mjjm(hi) is a zero-order generating functional. Using the Poisson summation 
formula Mjj(i)(hi) can be represented as 



x exp 



■*£<"?+»•£ 



M u{1) {hi 



ui{r x - r x+e 



n 



2vr 



+ 2ni>ff3 1 ^2r x m x + y^ihi 



(50) 



In addition to the above perturbation one has to extend the integration region over 
dual link angles uj\ to infinity. Clearly, all the corrections from this perturbation 
go down exponentially with (3 (it is obvious from the expression for the generating 
functional, Eq.(|5U|)). Calculating all the integrals in (|3U)l we find up to a constant 
(sum over all repeating indices is understood) 



M u(1) (h l ) = e$ htG u' h 'Y t e 



2 f3m x G xx /m x ,+Tr^ r i3hiDi(x')m x i 



(51) 



where we have introduced the link Green functions Gw and Di(x) (see Appendix 
A). G XjX > is the Green function of free massless scalar field in 3D. Up to the first 
exponent, the generating functional coincides with the dual partition function of 
the U(l) Villain model in the Coulomb gas representation. From here it follows: 
Only the configuration m x = for all x contributes to the asymptotic expansion 
and non-trivial monopole configurations are exponentially suppressed. We thus get 



Mum{hi) = exp 



- ^ hiGu'hi' 



(52) 
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Substitution of the last equation into (jHH) generates a weak-coupling expansion in 
inverse powers of f3. 

The easiest way to construct the corresponding expansion for the Wilson loop 
is the following. Let S^ y be some surface dual to the surface S xy which is bounded 
by the loop C and consists of links dual to plaquettes of the original lattice. Let b 
denote links from S^ y . Then, the Wilson loop in representation J is defined as 



W(C) = ( exp 



best} 



VP 



JuJh 



>■ 



(53) 



The last formula suggests that for the computation of the Wilson loop it is sufficient 
to make the shift hi —> hi + (i/ \/P)JJ2b^i,b in the expression for the generating 
functional (|5^). We find 



W(C) = exp 



'4/3 



Gbb' 



b,b'dSi y 



z u\i) n 



1 V- 1 A f 92 



exp 



- hiGwhy h — i=jhi 



Gi 



(54) 



It is now straightforward to calculate all connected pieces contributing to the Wilson 
loop at a given order of 1//3. For the first coefficients we get 



W(C) = exp 



where 



2P{C)= °w 

b,b'dSt, 



(55) 



(56) 



Using some properties of Gw described in the Appendix A it is easy to prove that 
this result coincides with the result given in for the U(l) model and for J = 1. 
The asymptotic behaviour of P(C) is given by Eq. (jHSJ)- 



3.2 577(2) LGT 

In this subsection we derive the low-temperature expansion for non-abelian models. 
Again, the general procedure is precisely the same as in 2D non-abelian spin mod- 
els, therefore for some technical details we refer the reader to paper [21]. First, we 
describe the general procedure for pure SU(2) LGT and then give a simple gener- 
alisation for all SU(N) models and include dynamical fermions. For simplicity we 
work in what follows with the Wilson action, i.e. H[Vi] = Re TrVJ. 

We want to expand the partition and correlation functions into an asymptotic 
series whose coefficients {B^g are calculated in some ensemble G with a Gaussian 
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measure. For the partition function we expect up to a constant factor an expression 
of the form 

oo . 

z = 1 + J2^ {Bk)G - (57) 

fc=l ^ 

Let us consider the standard parameterisation for the SU(2) link matrix 



Vi = exp[i<j k u} k { 



(5S 



where <r fc , = 1,2,3 are Pauli matrices and introduce 



W, 



1/2 



1/2 



where u> k (x) is a site angle defined as 



V x = exp [ta k u k (x)} = ( J[ V p ) C ( J] ^» ) C f • 



(59) 



(60) 



Here, the site x is dual to the cube c with vertices A and 5. us k (x) has the following 
expansion in powers of link angles 

u k (x)=^\x)+uW(x)+J*\x)+-- - . (61) 

On the dual lattice (we keep notations ()28|) for dual links) the first coefficient can 
be written down as 



uf\x) = uj k (k) + uj k (l 5 ) + u k (l 6 ) - w fc (Zi) - uj k (l 2 ) - Uk(h) 



(62) 



u% (x) for i > 1 can be computed by making repeated use of the Campbell-Baker- 
Hausdorff formula. Then, the partition function (|27j) can be exactly rewritten to 
the following form (see Appendix A in the e-print version of 



J SU(2) 



n 



sin 2 W, 



exp 



2/3 £ cos 



n 



w x 

sin W x 



\\ £ J JJrfa fc (x)exp 

x m(x)=—oo k 



-i a k {x)uo k (x) + 27Tzm(x)a(x) 



(63) 



where a(x) = (J2 k a \{ x )) 1 ^ 2 - For the derivation of this representation we have used 
the Poisson resummation formula. In order to perform the weak coupling expansion 
we make the substitution 



u k (l) - (2/?)" 1/2 ^(0 , a k (x) - (2(3) 1/2 Mx) 



(64) 



19 



and then expand the integrand of (|63|) in powers of fluctuations of the link fields. 
We introduce now the external sources h k {l) coupled to the link field u> k (l) and s k (x) 
coupled to the auxiliary field a k (x) and adopt the definitions 



u k (l) 







a k (x) 



d 



(65) 



dh k (l) ' ' ' ds k (x r 
With this convention we get the following expansion for the partition function IjfiHj) 



OO j 

1 + E jk B * (ft> 5 - 



fe=l 



M(h,s) 



(66) 



where the operators B k are defined through 



oo 

1 + E^ s *(ft.ft)=n 



X 



fc=l 



n 



9=1 



S (2/3)* 




u k >{x) 

, (2/3 )n/ 2 



■c fc w; 



2k 



(67) 



Here we have denoted 



a k 



2(fc+l) 



(2fc + 2)! 



> Cfe - E 



i 



n=0 



(2n + l)!(2Jfe-2n + l)! 



(6£ 



The first brackets on the rhs in Eq. (j67|) represent the expansion of the Jacobian. The 
first and the second brackets in the second line come from the expansion of the action 
and of the invariant measure, correspondingly. We have omitted the expansion of 
the term W x / sin W x since it does not contribute to the asymptotic expansion (due 
to the constraint W x = 0) but only to exponentially small corrections. As usual, one 
has to put h k = s k = after taking all the derivatives. Note, that in writing down 
these general expansions we do not distinguish between four BCC sub-lattices since 
the generating functional (see just below) has a unique form for all lattice. Thus, 
the product over x in Eq.()67|l goes over all sites of the dual lattice. The difference 
between BCC sub-lattices appears only in the general expressions for coefficients 
ou k l \x) for i > 1. 

The generating functional M(h, s) is given by 

p oo p oo i 

M(h,s) = / J^da fc (x) / JJdw fc (Z) exp[--u^(i) - iu k (l)[a k (x + e„) - a k (x)\] 



x.k 



l,k 



x exp 

m(x)=— oo 



2vri v / 2/3^ 



m\x )a\x 



x.k 



l.k 



(69) 
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where we extended the integration region for the link angles u) k (l) to the real axes. 
As in the abelian case, only the configuration with m x = for all x contributes to 
the asymptotic expansion, all others are exponentially suppressed. Calculating all 
Gaussian integrals in (joTJj) we come to 



M(h, s) = exp 



1 i 1 

-s k (x)G X}X >s k (x') + -s k {x)Di(x)h k (l) + -h k (l)G w h k (V) 



(70) 



From the last expression one can deduce the following simple rules 



Giv , (a k (x)a n (x')) 



>kn 



G xx , , -i(uj k (l)a n (x')) 



(71) 

The expansion the representation ()7()jl for the generating functional and the 
rules fl71j) are the main formulae of this section which allow to calculate the weak 
coupling expansion of both the free energy and any fixed-distance observable. The 
extension of this expansion for the Wilson loop is straightforward. Since all gauge 
invariant quantities depend only on the dual link variables but not on auxiliary fields 
one can make their direct expansion in powers of link angles. 



3.3 SU(N) LGT and dynamical fermions 

A generalisation for an arbitrary SU (N) model can be done again precisely like in 
2D spin models. For completeness, we repeat the arguments of [21] below. It is 
seen from the procedure outlined above that the low-temperature expansion arises 
only from the "vacuum" sector with m x = for all x. It follows from Eq. (J63j) that 
in this sector the SU(N) delta-function reduces to the Dirac delta-function so that 
the partition function becomes 



U dVi 



exp 



13 Re TlV i 



N 2 -l 



II I [ da k (x)exp[-ia k (x)u k (x)} . 



x k=l 



(72) 

After this simple observation the expansion itself is done precisely like for SU(2). 

To include fermions in the above expansion one can proceed in a standard way. 
Namely, we expand gauge matrices entering the fermionic matrix (J37)) as 



2<« 

U n (c) = exp[it k u k (c,n)] = } —[t k u k (c,n)] 

L — 4 ml 



m=Q 



(73) 



Here u k (c,n) = u) k {p) can be calculated from equations (|23|) and (|2^)l with the help 
of the Campbell-Baker-Hausdorf formula. We introduce now the following matrices 



1 3 



(74) 
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and 



■m 

F(c(x),c(x');U l [V l ]) = J2— E F^ m \c(x),c(x'))t kl ...t kn 



ml 

m=l k\...k m 



where 



I 2 

F {m) (c(x),c(x')) = -^[5 c+en!C ^ n (c)u kl (c,n)...u km (c,n) 



n=l 



+ (-l) m 6 c - entC ,Uc>k 1 (c',n)...uj km (c',n)] 
Then, the fermionic part of the partition function (J27|) is expanded as 



exp 



J^Trln^ [c{x),c{^)\Ui\yi 
/=i 

" s 

x exp 



V oo 

EE 

f=l 8=1 



IjDet^ 1 ] 



-Tr(GfF) s 



(75) 



(76) 



(77) 



With the substitution (jMj) and expressing the plaquette angles co k (c,n) in terms of 
w k (l) the last formula can be expanded further in powers of fluctuations of dual link 
variables. This procedure obviously does not change the generating functional there- 
fore all expectation values can be obtained again by differentiating the functional 

(HOD. 

In order to show how our expansion works in practice we have computed the 
simplest possible quantity, namely the first order coefficient in the expansion of 
the free energy. All details of these calculations for SU(N) model are given in the 
Appendix C. 



3.4 Some remarks on the low-temperature expansion 

We close this section with brief comments on general properties of the weak coupling 
expansion in non-abelian models. Consider the zero order partition function. As 
follows from Eqs. (j5U|) and (J72|) for the SU(N) group it has the form 



Y\du k (l) exp 



i,k 



I.k 



n 

x,k 



5\J k °\x) 



N 2 - 1 . (7£ 



The integrations are performed over the real axes though originally they are re- 
stricted to a compact domain, e.g. for 577(2) to (after change of variables flMJ)) 
Y^ fc u;f(Z) < tt 2 (3. From here it follows that the probability for K links to have large 
fluctuations like ui k (l) ~ 0(y/]3) is suppressed as 0( e - cori ^KTT 2 i3y This j s Q f course 



22 



only a rough estimate, nevertheless, qualitatively it is true and shows that large fluc- 
tuations are under good control since they are exponentially suppressed with /3, and 
this obviously remains true in the TL. As we have already stressed, this property is 
an essential achievement of the expansion in the plaquette representation. No such 
control exists for large fluctuations of the link gauge matrices in the large volume 
limit. 

The next observation which can be deduced examining the zero order partition 
function is that the low-temperature expansion is an expansion which starts from 
the abelian Bianchi identity for each component of the non-abelian field and goes to- 
wards restoration of the full non-abelian identity (compare with the strong coupling 
expansion, 12 ). 

It follows from the definition of the link Green function (|83|) that 

\G w \<4/3. (79) 

Therefore, with respect to the Gaussian measure in Eq.(|5U|) the fluctuations of the 
link variables are bounded by 

| (mu v ) G 1=1 \g w |< ~ , (80) 

due to the bound (J79|) . One sees from the last formula that the interaction be- 
tween links (=original plaquettes) strongly decreases with the distance. Taking the 
asymptotic expansion for the Green functions (j96|) we find 

| (^> G |=0(iT 3 ) . (81) 

This property justifies the low-temperature expansion in powers of fluctuations of 
plaquette variables, while there is no such justification for the expansion in terms of 
the original link variables fluctuations which are not bounded when L — ► oo. 

The zero order partition function (J7%|) is equivalent to the (iV 2 — l)-component 
free scalar field ak(x). In our expansion it enters as auxiliary field. As is seen from 
the generating functional (fTOj) it depends only on infrared finite Green functions. 
Link Green functions are infrared finite by construction (see Appendix A), while 
the Green function for the free scalar field is finite in any dimension D > 3. For 
D = 3 one has in the TL 

I Gw |< G = 0.5634 . (82) 

Nevertheless this shows that while correlations of dual links and correlations of links 
with auxiliary fields decrease with distance as is seen from ()71|) this is not the case 
for the correlations among auxiliary fields. It follows from f!71|) and the last bound 
that this correlation stays constant however large the distances between fields are. 
Just this property together with certain contributions from connectors creates an 
infrared problem in non-abelian models. To be precise the problem appears in the 
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following manner. Connectors generate contributions of the form Y^iec u k{l) and of 
high order to all u$(x), see Eas. (|innjl . (fTnHD - ffTT77j) . And though all u k (l) behave 
like 0(l/y/P) it is not obvious if it remains true for the sums of the plaquette angles 
along connectors. Since in the TL practically all connectors become infinitely long, 
this raises the question whether the property Y^iec^kil) ~ C(V\/7^) holds in the 
limit L — > oo? 

4 Conclusions 

In this article we proposed a plaquette formulation of non-abelian lattice gauge 
theories. Our approach to such formulations is summarised in section 2.1. We have 
also included dynamical fermions in our construction. The main formula of section 
2, Eq. (j27|) . gives a plaquette formulation for gauge models with dynamical fermions. 
As an application of our formulation we have developed a weak-coupling expansion 
which can be used for a perturbative evaluation of both the free energy and gauge 
invariant quantities like the Wilson loop. We believe that this work can be useful 
in at least two aspects. The first concerns the problem of the uniformity of the 
perturbative expansion in non-abelian models and is described in section 3. The 
second aspect concerns the perturbative expansion of lattice models with actions 
different from the Wilson action. Practically, all standard actions discussed in the 
literature have a very simple form in the plaquette formulation. The hardest part 
in the perturbative expansion is to treat contributions from the Jacobian. Since 
the Jacobian represents Bianchi constraints on the plaquette matrices and is the 
same whatever original action is taken it is sufficient to compute contributions from 
the Jacobian only one time and use the result for all actions. For example, the 
expressions for the coefficients C meas , Cj\ and Cj2 in (jlf)9|) given by (jlll|) . ()119|) 
and ()12()j) . respectively must be the same for all lattice actions. The expression for 
C ac may vary from action to action only. 

In our next paper j^Zj we derive an exact dual representation of non-abelian LGT 
starting from the plaquette formulation and study in details its low-temperature 
properties. In particular, we shall compute low-temperature asymptotics of the 
dual Boltzmann weight, derive its continuum limit and obtain an effective theory 
for the Wilson loop. 

In conclusion, we hope that the present investigation gives a certain background 
for an analytical study of gauge models in the low-temperature phase which is the 
only phase essential for the construction of the continuum limit. It can give a solid 
mathematical basis for conventional PT by proving (or disproving) its asymptotic- 
ity in the large volume limit. We also think that the present method can give 
reliable analytical tools for the investigation of infrared physics relevant for the non- 
perturbative phenomena like quark confinement, chiral symmetry breaking, etc. 
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A Properties of the link Green functions 



In this appendix we would like to describe the basic properties of the link Green 
functions Gw and Di(x') which appear to be the main building blocks of the low- 
temperature expansion in the plaquette formulation both in abelian and in non- 
abelian models. These functions appear also in the low-temperature expansion of 
the 2D spin models in the link representation. They are defined in precisely the same 
way as in 2D and have many properties analogous to properties of link functions in 
2D, see [21] for details. The functions Gw and Di(x') are defined as 4 

Gw = 25lj> — G X)X i — G x+en)X i + e nl + G X)X i + e n , + G x+enyX i , (83) 

Di(x') = G x>x i — G x+eniX > , (84) 

where the link / = (x, n) is defined by a point x and a positive direction n. G X)X > is 
a "standard" Green function on the periodic (or some other) lattice 

Gx,x> = jsYl 7(k) ' f( ~ k) = 3 - ^ C0S T kn ' ^ ° ■ (85) 

k n =0 ^ ' n=l 

The normalisation is such that Gu = 4/3. In the momentum space Gw reads 

Gw = 28„, - ± m (1 - e- fc ")(l - e-- fc "') • (86) 

Using this representation it is easy to prove the following "orthogonality" relations 
for the link functions 

Y^GibGw = 2G U > , ^2D b (x)G w = , ^ D ^) D bW) = 2G X , X , , (87) 

b b b 

where the sum over b runs over all links of 3D lattice. Let O be any closed surface. 
Then 

J2 Gw = , (88) 
i,v en 

where Gw = Gw if both link I and I' point in either positive or negative direction 
and Gw = —Gw if one (and only one) of the links points in negative direction. In 
particular, Eq. (j88|) holds for each cube of the original lattice. 

Let li, i — 1, 6 be six links attached to a given site x (notations are as in (J28j) ). 
One sees that Gw satisfies the following equation 

Gi t i> + Gi 2 i> + Gi s i' — Gm> — Gi 5 i> — Gi 6 i' = (89) 

4 Please, note that what we call link Green function does not correspond to link Green function 
of [2>| but rather coincides up to a Kroneker delta with plaquette Green function of Eq.(2.24). 
We call them link functions only because in 3D plaquettes are dual to links. 
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for any link I'. Di(x') obeys the lattice Laplace equation 



D h (x') + D h (x') + D h (x') - D u (x') - D h (x') - D h (x') = 26 x>x , . (90) 

The last three equations ensure the surface independence of the Wilson loop in the 
plaquette formulation, they allow to deform some given surface to any other one. 

Let S!L be the surface dual to the surface spanned by the Wilson loop. Consider 
the Wilson loop C with sizes R x T and lying in the x — y plane. Then dual links 
from the minimal dual surface S xy have the following coordinates I = (xi, 22, ^3; n 3 ), 
where x% is arbitrary and < Xi < R — 1, < X2 < T — 1. We then have for the 
quantity P{C) entering the weak-coupling expansion of the Wilson loop (J55|) 

£ Giv = 2P(C) , (91) 

where 

nfn\ 1 2 - cospi - cosp 2 1 - cospii? 1 - cosp 2 T 2tt 

L 6 t—* 1 j(k) 1 — cospi 1 — cosp 2 L 

In the TL and for R, T —>■ 00 one finds for 3-D 

P(C) ^ -(RlnT + T\nR) , (93) 

7T 

what coincides with [26 . 

Finally, we study the asymptotic properties of the link function Gu>. Let I = 
(0, X2, £3; rii) and I = (R, x%, £3; n\). From the definition (J%3~j) it is easy to obtain 

G w = 2D{R) -D(R + 1) -D(R-l) , (94) 

where 

1 x — v 1 cos j kv-> -Rj'n 

m) = j; E — ^ — • < 95 > 



Since when R — > 00 



we come to the following asymptotic behaviour 



„ const . , 

Gw * • (96) 
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B Relation between site and link angles 

Define the site angle uj k {x) as follows 

V x = exp [it k u k {x)\ = \JJjtM)J e *Mc) (jjjt k euk)\ e -it k e k (C) ? ^ 

where V x G SU(N), t k are the generators of the SU(N) group and 

e it k 8k(C) _ e it k d k (b) _ ^gg^ 

bee 

oo k (x) has an expansion in powers of the link angles 6 k (l) which can be obtained by 
making repeated use of the Campbell-Baker-Hausdorf formula. The expansion is of 
the form 

io k (x) = J°\x) + w£\x, C) + J k \x, C) + ■ ■ ■ , (99) 

where we have indicated explicitly the dependence of the coefficients on the connec- 
tor C. Let us denote 

4 0) (^) = 5>w • ( 10 °) 

bee 

Then the first three terms in the expansion (|9*9*j) have the following form 

^\x) = j^e k {k) , (ioi) 
i=i 

J k 1 \x,C)=J k 1 \x)+J k 1) (C) 1 (102) 

^\x,C)=J k 2 \x)+J h 2 \C), (103) 

where we have separated the connector contributions explicitly, i.e. oo k \x) depend 
only on six links U which belong to a given site x, while u k {C) depend also on links 
b G C. All these coefficients can be written down as 

6 

<^k\ X ) = -tkmn e m(h)0 n (lj) , (104) 

i<j 



1 6 

6 

3 



3 

ij=l 



+\ Yl [20n(ii)0fc(/i)0n(/ P )-0«(ii)0n(/i)0fc(ip)-0fc(/i)^(/i)0n(ip)] , (106) 



i<j<p 
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+2£ 



i=4 



3 6 

i=l i=4 
6 

3<i<j 

Y, e k(b)e n {b)e n (u) + 2 e n {b)9 k {b')e n {k)-9^\c)e^{C)e k {k) <io7) 

6GC 6<fe'£C 



C Exact expression for the first coefficient of the 
free energy expansion of SU(N) model 

Here we present the computation of the first order coefficient of the free energy 
expansion of the SU(N) model 

Nf 

F=^lnZ = P-^^lnP-y^^L + ^--C 1 + 0(p- 2 ) . (108) 
3iVL 3 3N 7^ 3Ar ^ 3 3Nf3L 3 Kf ' V ; 

There are five contributions to C\ 

C\ = Cac + Cmeas + Cjl + Cj2 + Cferm ■ (109) 

The contribution from the action is given by 

_J_ (iV 2 -l)(2iV 2 -3) 1 ~ 2 (jV 2 -l)(2iV 2 -3) 

3NL* ac 96iV 2 3L 3 ^ H ~ 54iV 2 1 Uj 

and the contribution from the measure by 

C meaS = ~ ^Tg ^ G " = 18 ' 



3iVjL 3 — 24 3L 3 



As is seen from the first line in the Eq. (l67|) there are two contributions from the 
expansion of the Jacobian. They are given by the following expectation values 

1 On = ^H(~h £ U(x)u,W{x)) , (112) 



3iVL 3 J 3NL 3 2 

x.k 



1 „ 1,1 



3iVL 3 3iVL 3 



<4> (fe 

\ \ x,k 



2\ 



C * = ^(-7) ( I > 1 ) . (113) 



To compute these expectation values we have to use the representation of site angles 
cjfr (x) in terms of link angles given in Appendix B, where 6^(1) coincide with 
Uk(l) up to a sign. Therefore, we can define 

wfe(0 = vMi) , m = ±! ■ ( 114 ) 
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This sign rji can be easily established for any link from the explicit expressions for 
the Bianchi identities given in Eqs.([29 p -()36| ) . Let k, i = 1,...,6 denote six links 
attached to a site x(c) from the c-th BCC sub-lattice. Let C be a set of links b 
which form a connector for a given site x(c). Introduce now the following quantities 

1 6 

Mi(a?(c);ZiZaZ 3 ) = [ 5 kh 8 kh 5 hh ~ khkh^h] 

ij=l 

2 6 

+ 3 E r lk%'% [ 2S Uh S ijiiKh ~ ^khhjlAph ~ Kh 5 ijh 5 lph] > ( 115 ) 

i<j<p 



3 6 

i=i j=4 bee 

6 

+ 2 Yl Wj^Vb Pw/iifa^Ns -^fa^Ijti^Ms] 

3<i<j feec 



i=4 L&eC 



E 7 & < WwA'a + 2 VkVbVb'SbhSyiAih - E VkVbVb'5 b i 2 5b'i 3 5i lh 



b<b'ec 



b,b'ac 



B(x(c),x'(c);hl 2 ) = B(x(c);l 1 )B(x'(c');l 2 ) 



(116) 
(117) 



B(x(c);h)= ^m.mAh + ^Vb^mM • (us) 

_i<j beC i=4 

With the help of these quantities, and after taking all the derivatives and calcu- 
lating the sums over the group indices we present the result in the form 

N 2 -l 



3NL 3 



a 



.a 



2N 



[A 1 + A 2 ] , 



3iVL 3 



32N 



where 



4 = E £ M *( x ( c MW 

C=l 1(c) W 2 «3 

[(iV 2 - l)D {l (a:(c))G /a/ 3 + D, a (a;(c))G Ws + D h {x{c))G hh ] . 



(119) 
(120) 



(121) 



i 4 

^ i = 3^E E E fl ^ c )'^i^)[ G wAW c H^ + (122) 



c,c'=l x(c),x'(c') (i^2 



29 



G ljjl2 D v ,Xx{c))D h {x'^)) - G h , l2 D lf/ (x(c))D h (x'(c')) - G ljA D h (x(c))D h (x' '(c'))] , 

3 J J 3 



® 2 = 3L 3 " ^ ^ ^-B(a ; ( c )7a; / (c / );^i^)G ? :!; ( c ) i:r '( c ') Gi 1 ^Gi hh - G^^G^^ 



c,c'=lx(c),x>(c') hh 



(123) 

Finally, we compute the contribution of fermions. As follows from (|75|) - (j77jl two 
l//3-contributions are given by 

Cferm = ~ EE^^^W) + ^ EE^WHIpH^)) ■ (124) 
/ P / PP' 



We thus obtain 



1 



T 2 - W / 



3NL 3 



C 



ferm 



N - 1 
48L 3 



5> 

/=! 



/ 



-Ew E 4Ew>EE 



(/,/')Gp pp' 



iep /'ep' 



where we have denoted 

H f (p) = Tr s G/(c + e n , c)f n (c) + G/(c, c + e n )£ n (c) 



(125) 



(126) 



(127) 



H f ipp') = Tr s [ G f {d + e n ,, c)Uc)G f {c + e n , c')^(c') 
+ G f (c', c + e n )£ n (c)G f (c, d + e n >)£ n ,(d) 
- G f (c', c)£, n (c)G f (c + e n , d + e n >)£ n ,(c') 

G f (d + e nl ,c + e n )Uc)G f (c,d)^(d)] (128) 

and Tr s means the trace over spinor indices. 
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